Non-equilibrium growth in a restricted-curvature model 
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The static and dynamic roughenings of a growing crystalline facet is studied where the growth 
mechanism is controlled by a restricted-curvature (RC) geometry. A continuum equation, in anal- 
ogy with the Kardar-Parisi-Zhang (KPZ) equation is considered for the purpose. It is shown here 
that although the growth process begins with a RC geometry, a structural phase transition occurs 
from the restricted-curvature phase to the KPZ phase. An estimation of the corresponding critical 
temperature is given here. Calculations on the static phase transition give results along the same 
line as existing predictions, apart from minor numerical adjustments. 



PACS number(s): 05.40.+j, 68.35.Rh, 64.40.Ht 

Growth of interfaces in a strongly temperature con- 
trolled regime has been the subject of a considerable por- 
tion of recent technological developments [1] . The inter- 
est in this field has mainly been generated by the micro- 
scopic roughness that originates as a result of competition 
among different effects, such as surface tension, thermal 
diffusion and different noise factors coming into play dur- 
ing the growth process. Even in equilibrium, a crystalline 
facet "remains practically flat until a transition temper- 
ature is reached, at which the roughness of the surface 
increases very rapidly" [2]. Numerous theoretical mod- 
els, starting from Burton, etal [2] have been proposed to 
account for a detailed analysis of the height fluctuations 
during the growth process on both sides of the roughen- 
ing temperature Tr [3-5] . They found a change in mobil- 
ity from activated growth in the nucleated phase (char- 
acterized by low temperatures) to nonactivated growth 
at large temperatures. The pinning-depinning growth 
model of Chui and Weeks [3] has later been numerically 
extended to a polynuclear growth model by Sarloos and 
Gilmer [6] confirming a growth by island formation which 
is continuously destroyed by any small chemical poten- 
tial favoring the growth. Theoretical forays in this front 
have continued with detailed predictions on the equilib- 
rium roughening transition [7,8] which have also found 
experimental justifications in [9]. 

Starting from the early theoretical efforts [2,3] to 
the present day developments [1,10], numerous discrete 
[11,12] and continuum [13,14] models have been proposed 
to study both equilibrium and nonequilibrium surface 
properties. However all these different models seem to 
fall within either of the KPZ or the Lai-Das Sarma uni- 
versality class, barring a few exceptions [11,15]. The 
atomistic growth process of the latter type deals with sur- 
faces grown by molecular beam epitaxy (MBE) method, 
whose distinguishing feature is that growth occurs under 
surface diffusion conditions, with the deposited atoms re- 
laxing to the nearby kinks. The essential idea employed 
was to modify the relaxation mechanism as a locally sur- 



face minimizing curvature, instead of surface area. To 
linear order, this mechanism was supposed to mimic the 
growth dynamics of crystalline surfaces [12,14]. 

Alternative efforts in this front have mainly centered 
around the development of theoretical models whose dy- 
namics can be mapped to cither of these two main uni- 
versality classes. A classic example is an equilibrium 
restricted- curvature (RC) model studied by Kim and 
Das Sarma [16]. The corresponding growth rule describes 
a growth restricted on the local curvature, | V 2 h |< N, 
and is obeyed at both the growing site and its nearest 
neighbors where N is any fixed positive integer. We 
adopt the nonequilibrium class of growth proposed by 
Kim and Das Sarma as the starting point of our study 
of the dynamics of the nonequilibrium growth in a MBE 
process. 

Considering two dimensional growth pertaining to a 
lattice structure where h(r, t) is the height of the interface 
at time t at position r, the equation goes like 



r)^ = -lV 4 Hr,t) + ±\V 2 h(r,t) f 

sin[—(r, t)}+ F + R(r, t) 



(1) 



where 77 is the inverse mobility which fixes the time 
scale, 7 is the surface tension, a is the lattice constant 
and V is the strength of the pinning potential. F is a 
steady driving force with R the white noise defined as 



< R(f,t)R(r , ,t') >= 2D5 2 (f- r')5{t - t') 



(2) 



The biharmonic term on the right side of the above 
equation gives the basic relaxation mechanism in opera- 
tion. The second term is the most important term for the 
restricted-curvature (RC) dynamics and has been incor- 
porated in analogy with the KPZ equation [13]. Just as 
- 2 

the lateral term | Vh \ gives the nonlinearity in a KPZ 
growth, the curvature dependence in our model is pro- 
posed to produce a leading order nonlinearity of the form 
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(V 2 h) 2 . From phenomenological considerations, since 
the discretisation scheme prohibits a high curvature, the 
constant A is negative. The third term comes along due 
to the lattice structure preferring integral multiples of h, 
which means that h is measured in units of a. F is the 
steady driving force required to dcpin a pinned interface 
while all the microscopic fluctuations in the growth pro- 
cess are assimilated in the noise term R. It can be shown 
easily that the second and third terms are obtained from 
a variant of the sine-Gordon Hamiltonian 

E[h(f,t)]=J J <Pr[l(V 2 hf-Vcos{ 2 -?-h)] (3) 

We now start pursuing the most important goal re- 
garding the dynamics of the growing interface - whether 
a non-equilibrium phase transition exists or not. In other 
words, whether we can specify a critical temperature 
Tr below which the surface is flat and above which the 
surface starts showing a dynamic roughening. This ob- 
viously brings into question the interplay of strengths 
between the non- linear term (V 2 /i) 2 and the sine- Gor- 
don potential. Three limits of the above eqn.(l) are well 
known: 

(i) A = V = F = 0; characterizing a stationery interface 
[11,12,16]. This sort of growth has found experimental 
justification in the works of Yang, etal and Jeffiies, etal 
[19]. 

(ii) A = F = 0, V ^ 0; characterizing the growth dynam- 
ics of crystalline tensionless surfaces [20] . This model in 
equilibrium depicts a roughening transition to the high 
temperature regime of the sine-Gordon model and is ex- 
pected to modelise the vacuum vapor deposition dynam- 
ics by MBE growth [21]. 

(iii) V = F = 0, A^0;a very special case of the 
numerical simulation predicting a " local model" for den- 
dritic growth with analytics failing to define a stable, 
non-trivial fixed point [22] . This issue is still under much 
debate and really needs a deeper understanding to have 
any final say in the matter. 

However in both the first and second situations, de- 
tailed numerical analysis have predicted a transition from 
the conserved MBE growth to the Edwards- Wilkinson 
[23] type phase characterized by the generation of a V 2 fe 
term [20,24]. The generation of this so-called "surface 
tension" term in the dynamics although surely being a 
fall-out of the renormalization of the sine-Gordon poten- 
tial in the latter model, actually demands a greater at- 
tention. Combining this idea with the proposition put 
forward by Kim and Das Sarma [16], we therefore pose 
the most general situation concerning a surface growing 
under MBE and try to ascertain the associated rough- 
ening process throughout the whole temperature range, 
with the growth essentially occuring under a surface cur- 
vature constraint. In the following analysis, we employ 
standard dynamic renormalization techniques to probe 



the dynamics of our proposed Langevin equation (1) [7,8] 
and later on follow the line of Chui and Weeks [3,7] in 
taking account of the static renormalization of the model 
Hamiltonian shown in eqn.(3). 

As usual, we start by first integrating over the mo- 
mentum shell A(l — dl) <| k |< A perturbatively in 
A and V . Thereafter we rescale back the variables in 
the form k — ► k' = (1 + dl)k, h — > h! = h and 
t — ► t' = (1 — Adl)t. The various coefficients follow the 
rescaling 77 — > 77' = 77, 7 — > 7' = 7, A — > A' = A, V — > 
V = (1 + Adl)V, F -> F' = (1 + Adl)F. We set up 
the perturbative scheme to rewrite eqn.(2) in a comoving 
frame moving with velocity Fj r\ as 

V^h = - 7 V 4 /i + $(/») + R (4) 

where = -^sin^ {h+^ t)] + f (V 2 h) 2 . There- 
after going exactly by the analysis of Nozieres and Gallet 
and Rost and Spohn [7,8] and including the corrections 
in [25], we finally arrive at an expression for the renor- 
malized mode coupling term, 

,so = _^ dl f_J d , rl _l_ Mmf _ n) 

God ?-? \,t~t') x e -[^(*-*') + !^(I^M-*')] 

(r, - r'j)} cos[~(t f )] + [1 - ^-[ di h(r,t)} 2 

x(n-rif} Sl nl~(t-t>)} (5) 

where A ~ 1/a is a suitably chosen upper cut-off with 
the Green's function G(x, t) given by 

and 

<j>(p,x) = £ f[l-J o{ k~p)]e-^ 2 ^ (7) 
with p = Ap and x = 7 ^~ 2 * - . 

Turning now to the above eqn.(5), we find that starting 
with a structurally non-linear term in the Langevin equa- 
tion, the dynamic renormalization has initiated a mode 
coupling structure with quite a different composition. 
The absence of the nonlinear (V 2 /i) 2 term in $ SG implies 
that once starting with a restricted curvature model of 
growth, The lattice structure partakes a dynamics where 
in a finite time after the start, A is renormalized to zero 
and thereafter the KPZ type [13] nonlinearity takes over. 
Thus a competition ensues between the alternate pinning 
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and depinning forces offered by the (Vft.) 2 and (\7 2 h) 2 
nonlincarities. Also the production of the surface ten- 
sion term V 2 /i ensures a dynamic phase transition from 
the restricted-curvature regime to the KPZ regime. From 
an analysis of the following renormalization flows, we ar- 
rive at an expression for the temperature at which the 
transition occurs: 



(8) 



A (A,) (n;ft)= [ — [ dp p 3 J Q (p) x sin[ 
Jo x Jo 



2ir Kxp 2 
a 7 



X J d P e^ {rp '^^P Xp4) X e -[^ 2 +2n0(P,*)] ( 18 ) 



*(k\, ^ f°° dx f°° - t ,^ < 2t < 
A y >(n;K)= / — / dp p Jo(p) x svn[ — 

Jo x Jo a 



2ir Kxp 
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x j, e iM^ xe -^ (19) 



^ = ^! nj4 (7') (n;K) U2 

dl -fa 4 V ' ' 



dj it 4 
~dl ~ 67a 4 



dry 



8tt 4 r? 

7a 4 7 



iA(T>(n;/s) C/ 2 



nj4 ( ">(n;/s) U 2 



d\' _ 8tt 5 
d7 7a 5 



nA( A ')(n; K ) U 2 



(9) 



(10) 



(11) 



(12) 



(13) 



= 4^ + a - 14 ^ w (»; «) ^ 2 (14) 

at 47rr77 7a d 



dD 1 _ \ 2 D 8tt 4 D ... . rr2 
at 87r 7^ 7a 4 7 

where U = X, AT = k = ^ and n = with 
7' representing the coefficient of the renormalized surface 
term [26] and A' denoting the coefficient corresponding 
to the KPZ nonlinearity generated on account of renor- 
malization. A^> (n; k) stand as shorthand representation 
for the integrals, where i = 77, A', A, k, etc. and whose 
detailed functional forms are given below: 



A^'\n-K) 



dx 



i~ ~3 t ( ~\ r 2?r Kxp 2 

dp p Jo{p) x cos[ J 

' 07 



f°° dr f°° 

A^(n;K)= / — / dpp 5 Mp) 
Jo x Jo 



AW(n-,K) 



x cos 



2tt Kxp 

a 7 



/ 



dp e 



f(*p- 



x e 



-[xp 2 +2n<j>(p,x)] 



(20) 



As already argued, the eqn.(13) giving the flow for the 
renormalized KPZ coefficient generates the lattice poten- 
tial of the driven interface after the structural phase tran- 
sition has occured. At some temperature Tr = -^-,in 
the units measured, corresponding to the transition point 
n = 4, the roughness dynamics driven by the KPZ force 
takes over. From this point the whole dynamics follows in 
the line predicted by Nozieres and Gallet [7] and Rost and 
Spohn [8], with the biharmonic and structural nonlinear 
terms muffled by the surface tension term and the KPZ 
nonlinearity respectively. This KPZ regime continues un- 
til the growing interface reaches the next crystalline facet 
whereon the RC dynamics again comes into play and 
the whole mechanism goes on repeating itself through- 
out the process of nonequilibrium growth. However in 
the static case, starting with the standard sine-Gordon 
Hamiltonian (eqn.(3)) and again taking clues from [7,8], 
we arrive at the following Kosterlitz-Thouless type sec- 
ond order equations [5,10] in terms of the reduced vari- 
ables y = M and 1 = 



/*el<-^V, xe -«,.» (16) 



A r A~~T(~\ ^Kxp 2 

dx I dp p Jo(p) x cosy 

Jo « 7 



where 



and 



B{n) = f C 
Jo 



dp p 6 J (p) e 



-2nh(p) 



HP) 



dk 
F 



[1 - Mkp] 



(21) 



(22) 



(23) 



(24) 



/ 



Ipxp ) 



with n = The roughening transition occurs at 



x e [xp +in<t>{p,x)\ ( 17 ) tne fixed point y = q x = 1 and we recover a n tne pre . 

dictions of Nozieres, [10] although with slightly different 
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coefficients. The point to be noted is that although the 
transition temperature is twice here compared to that 
of the standard Kosterlitz-Thouless (KT) case, both the 
dynamic and static phase transitions occur at the same 
critical temperature Tr = somewhat alike to the 

NG and KT models. 

In conclusion, we have studied interface growth models 
both in the equilibrium and nonequilibrium cases under 
a constraint of the growth occuring under a curvature 
restriction. The dynamic model gives a phase transition 
from the RC growth to the KPZ type at a temperature 
we have calculated to be Experimental data on 

the roughening transition of solid 4 He in contact with 
the superfluid 3 He [9] seems to show a crossover to the 
non-equilibrium regime in the superfluid phase. We sus- 
pect that in a scenario analogous to the MBE growth 
model proposed, the crossover situation is actually mim- 
icked by the dynamic phase transition found in the model 
discussed and the larger value of the critical phase transi- 
tion temperature reported here might just be the expla- 
nation for the weaker coupling observed around 0.1°K 
in 3 He. However, although the dynamic mechanism pro- 
posed here differs largely from the other existing models, 
thereby defining a new universality class, the behavior 
at equilibrium defined by the corresponding static model 
is actually a replica of the Kosterlitz-Thouless criticality, 
the only variation being in an augmented value of the 
transition temperature. 
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